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Abstract 

The LHC has the potential not only to discover supersymmetry (SUSY), 
but also to permit fairly precise measurements of at least a portion of the 
sparticle spectrum. Proposed mass reconstruction methods rely upon either 
inverting invariant mass endpoint expressions or upon solving systems of mass- 
shell equations. These methodologies suffer from the weakness that one certain 
specific sparticle decay chain is assumed to account for all the events in the 
sample. Taking two examples of techniques utilizing mass-shell equations, it 
is found that also applying wedgebox and hidden threshold (HT) techniques 
allows for the isolation of a purer event sample, thus avoiding errors, possibly 
catastrophic, due to mistaken assumptions about the decay chains involved and 
simultaneously improving accuracy. What is innovative is using endpoint mea- 
surements (with wedgebox and HT techniques) to obtain a more homogeneous, 
well-understood sample set rather than just using said endpoints to constrain 
the values of the masses (here found by the mass-shell technique). The fusion 
of different established techniques in this manner represents a highly profitable 
option for LHC experimentalists who will soon have data to analyze. 
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1 Introduction 



LHC experimentalists will soon determine whether or not Supersymmetry (SUSY) 
is a TeV-scale phenomenon: if so, colored sparticles will probably be the first to be 
discerned, possibly soon followed by neutralinos, charginos, and sleptons if favorable 
decay channels are open, though measuring the masses of these latter colorless sparti- 
cles with percent precision will be challenging |1] |2] . The reason for this is that every 
R-parity-conserving^| SUSY event produces at least two invisible particles (the light- 
est SUSY particles, or LSPs) which carry away significant missing energy and make it 
impossible to reconstruct mass peaks. Therefore, many SUSY mass extraction tech- 
niques depend on precise measurement of invariant mass distribution endpoints. For 
a sparticle decaying into an LSP and a Standard Model (SM) fermion pair, either via 
a three-body decay or sequential two-body decays, it is straight-forward to see how 
the endpoint of the di-fermion invariant mass distribution yields the mass difference 
between the decaying sparticle and the LSP (perhaps modified by the on-mass-shell 
intermediate for two-body decays) [U [5]. Studies attempting to fully reconstruct 
the actual sparticle masses from invariant mass endpoint information rely on specific 
longer decay chains, typically q — > X2I ~* — ► Xi^ + ^~Q [HI El EE QUI [EE] 

each event would have two g's (for instance) produced. It is then theoretically doable 
to construct enough invariant mass distributions to determine the sparticle masses; 
however, in practice endpoint measurement may be complicated by low event rates, 
fitting criteria, unaccounted-for (in the simulation) higher-order and radiative [T2] 
effects, and (in particular) backgrounds. 

Ideas on how to measure SUSY masses without relying on distribution endpoints 
have also been put forward. The work of Nojiri et al. [T3l [H], for example, uses 
mass-shell relations in a sufficiently long SUSY decay chain, e.g. g — > b b — > x^bb ~^ 
l^Pbb — > l^Pbbxi', ^ some of the masses are already known in this chain the others 
can in principle be found be solving mass-shell relations for a small sample of events^, 
which may in fact lie far from the endpoint. Another method, due to Cheng et 
al. [T5], starts from very similar looking mass-shell relations, but instead of assuming 
some masses and solving for the others, scans the whole mass-space for points where 
these relations are most likely to be satisfied. Both methods, hereafter designated 
as 'Mass Shell Techniques' (MSTsjf], seem quite effective in obtaining percent-level 
determination of the sought-after SUSY masses, at least at the parameter points 
considered in those works. 

The accuracy of both these MSTs hinges on one critical assumption: the decay 
topology of choice has been isolated. In the actual LHC data, the decay topology 
would have to be inferred, if this is at all possible, before proceeding; MSTs would 
thus appear to be excellent roads to SUSY mass reconstruction which, however, begin 
only at a point half-way to the destination. 

The present work focuses on the first half of this road; i.e., isolating a desired 

3 If R-parity is not conserved, then it may be possible to fully reconstruct events. See [3J for 
further details. 

4 In practice, many events are still required. 

5 Refs. [131 HH HH] use the name 'mass relation method' to refer to their technique. 



decay topology, and on how this affects a subsequent MST analysis. As a first foray 
into this potentially quite thorny task, let's restrict ourselves to just the specific 
topologies considered in p2] and [15] involving a pair of neutralinos XiX'j (hj — 
2, 3, 4) that subsequently decay to leptons (electrons and muons) via on-shell sleptons. 
This situation is amenable to a wedgebox analysis [T71 [18] which is based upon a 
scatter plot of the di-electron mass M e + e - versus the dimuon mass M M + M -. A key 
benefit of this technique is that it allows (at least partial) separation of individual 
events according to the specific (i, j)-pair whose production gave rise to them. Given 
sufficient statistics, events from each such decay- type fall in distinct, easily-recognized 
zones of the wedgebox plot. The overall topology of the resulting wedgebox plot 
then tests for the significant presence of the various possible decay channels - 
which may for instance signal the meaningful presence of a decay channel erroneously 
assumed to have been insignificant as the basis for a MST analysis. Events can be 
selected from a specific zone of the wedgebox plot, preferably a zone populated by 
only one decay channel. This acts to maximize sample homogeneity and assure the 
basic MST assumption is satisfied. 

Although the wedgebox technique relies of locating the endpoints of invariant 
mass distributions — just like the studies 0, El El El QUI EH] mentioned earlier, the 
information sought is radically different: the wedgebox analysis is tailor-made for 
event sample sets comprised of assorted produced sparticle pairs and multiple spar- 
ticle decay chains. The observed endpoints serve to delineate the zones and allow 
for selection of purer subset(s) from an overall sample set. (Using this endpoint 
information to determine a set of cuts is a far more rational course than that of arbi- 
trarily choosing some numerical cut-off values to purify the data sample.) Virtually 
all previous studies presume such purification has already been accomplished either 
by an unspecified set of cuts or a fortuitous choice of SUSY input parameters. The 
wedgebox technique illustrates a concrete method of how to deal with the more gen- 
eral case, and the consequences which result and should not be ignored in a coupled 
analysis aiming to extract the sparticle masses. 

While the wedgebox plot structure does strongly indicate neutralino pair decaj|§, 
the topological structure seen (such as a box, a wedge, a double-wedge, or a box 
within a wedge for instance) does not differentiate between neutralino decay chains 
with either on- or off-shell intermediate particles, though the numerical values of the 
endpoints may show notable sensitivity to this [I]. While the di-lepton distribution 
shapes for on- and off-shell decays may be theoretically distinguishable 0, [17] , even 
if this is true the distinguishability may be lost after factoring in the effects of cuts, 
backgrounds, etcQ If the neutralino decay is via an off-shell intermediate, then use 
of an MST fails simply because some of the mass-shell conditions fail. The recently 
discovered Hidden Threshold (HT) technique, however, can test the on- or off-shell 
nature of a sample of decays [19]. Here the idea is to look for correlations among 

6 The only other possibility in the MSSM involves a chargino decaying through a sneutrino, 
X2 —> —> l + l~xf> since the kinematics are the same. The wedgebox plot must be considered 
effective in separating these decay topologies as well. 

7 This same point was made in [5]. See Discussion Point # 2 and Section 5.2 therein. This will 
be examined more closely in Seen. 



various invariant mass combinations in the four-lepton endstate - a genuine off- shell 
decay fulfills the correlations while an on-shell decay does not. 

This paper will show that MSTs, by construction unrelated to invariant mass 
endpoints, can nevertheless be improved by information contained in these endpoints 
- specifically via the wedgebox and HT techniques^]. The paper is organized as 
follows: Sections [2] and [3] concentrate on the 'N-MST' method of Nojiri et al. and 
the 'C-MST' method of Cheng et al., respectively; Section [5] then offers conclusions 
and some additional discussion on SUSY mass spectral analyses at the LHC. 



2 The N-MST Method of Nojiri et al. 

For the N-MST method, the focus will be on the decay of a heavy MSSM pseudoscalar 
Higgs boson as considered in [13]. The specific decay chain is 

pp - a° - x°x° - 5 VSV - Zi^iV^xiXi (i) 

where the Higgs boson decays to neutralinos (i, j = 2, 3, 4) via on-shell sleptons of the 
electron- or muon- variety (see Fig. [T]). Assuming the final leptons' four- momenta 




Figure 1: Feynman diagram for the decay Here i,j = 2,3,4 while the label 1 or 2 on 
leptons stands for either e or fi. 

are known while the LSP's escape detection, ([1]) implies six mass-shell constraints on 

8 Though our strategy is quite different from the 'hybrid' method of [16) which couples an MST 
with values of endpoints. 



the eight unknown components of LSP four-momenta (p M and 
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Nojiri et al. also posit two overall momentum conservation constraints 
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along directions (x , y) transverse to the beam (z), though this would appear to be 
contingent on the Higgs boson having no transverse momentum. If all the masses 
m~o, m,~, m,~, m~o, and m~o are known in advance, one can solve the eight equa- 
tions (J2J)- (J9j) for the eight unknowns and reconstruct the Higgs boson mass via 



from just one Higgs boson event of the type ([I]). However, even in this idealized sce- 
nario which does not include detector resolution, particle widths, backgrounds, etc., 
there are two major caveats which prevent this procedure from being so straightfor- 
ward: 

• There is a 4-fold ambiguity in assigning labels Zi^, l[ 2 to the leptons; this forms 
a combinatoric background. 

• Combining ©-Q leads to a quartic equation with 0, 2 or 4 solutions for the 
unknown momenta. 

So what one must do in practice is collect a number of events and deduce the correct 
value of ttla from the maximum of the resulting distribution. In [13J a 10 3 event 
sample (with no backgrounds) thus yielded a percent-level determination of the Higgs 
boson mass. 

2.1 Addition of the Wedgebox Technique 
2.1.1 Box Topology 

As shown in [13] the programme sketched in the previous section works fairly well 
at Snowmass Benchmark SPSla [20j where the dominant Higgs boson decays are 
via ([1]) with i = j = 2. Therein, 1000 events of the type A — > X2X2 ~^ MHXiXi 
were generated using the HERWIG6.4 [2T] event generator and passed through the 



m A = {Ph + Ph + Pl[ +Pl' 2 +P + P'f 



(10) 



detector simulator package ATLFAST [22]. The only cuts required were that all 
four isolated leptons have r\ < 2.5, with two of the leptons also having p T > 20 GeV 
while the other two have px > 10 GeV. Same-flavor events such as e + e~e + e~ or 
were also included if one of the two possible pairings of OS leptons in 
such configurations gave a di-lepton invariant mass beyond the ~ 78 GeV kinematic 
endpoint (implying the other possible pairing is the correct one). This analysis 
yielded the correct Higgs boson mass with a resolution of only 6 GeV. 

The following set of cuts are herein adopted in a effort to reproduce these results 
using IS A JET [23] in place of HERWIG: 

1. Leptons must be hard(p T > 10, 8 GeV for e ± , /x ± , respectively; |?/| < 2.4), 
and isolated (no tracks of other charged particles in a r = 0.3 rad cone around 
the lepton, and less than 3 GeV of energy deposited into the electromagnetic 
calorimeter for 0.05 rad < r < 0.3 rad around the lepton). 

2. There must be missing energy in the range: 20 GeV < E)t < 130 GeV . 

3. Nojetfl are present with a reconstructed energy Ej et greater than 50 GeV. 



A sufficient sample of pp — > A events is collected to represent an integrated LHC 
luminosity of 300 fb~ . With this stricter set of cuts it should be somewhat harder 
to extract the Higgs boson mass. 

Fig. [2^ shows the wedgebox plot at SPSla for e + e~/z + /i~ events. A 'simple box' 
topology consistent with the expected X2X2 origin of lepton pairs is clear. Moreover, 
the number of flavor-balanced events (e + e~e + e~ + + e + e~fi + [i~) exceeds 

the number of flavor-unbalanced events (e + e~e ± /i T + / u + /i _ / u ± e =F ) at this parameter 
point [28]; this indicates that the events come primarily from a Higgs boson deca}0 
with decay topology ([T]). Though the final number of events passing cuts is somewhat 
small (only 30 or so compared to the 1000 which [13] assumes, mainly due to the 
inclusion of the leptonic branching ratios) the number of events in the distribution 
of solutions is quite a bit larger: recall the bulleted caveats earlier which potentially 
can yield a 4 x 4 = 16-fold multiplicity factor. Though only a factor of 4 or so is 
observed, nonetheless a fairly clear peak in the distribution emerges (see Fig. [2b) at 
the correct valud^ of rriA = 400 ± 5 GeV. 



2.1.2 Wedge Topology 

In principle, this method should work for the case i 7^ j as well; to test this we 
consider the following MSSM parameter point: 

9 Jets arc defined by a cone algorithm with r = 0.4 and must have \r]j\ < 2.4. 

10 Note that though in a simulation one can of course only choose to generate Higgs boson decay 
events, experimentalists lack this freedom. 

11 Also note here that the momentum-conservation constraints ([9|) do not appear to be generally 
true; i.e., the parent Higgs boson is often generated with significant transverse momentum in the 
range < pt ~ 100 GeV according to this analysis. Surprisingly, this does not seem to affect the 
result. 
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Figure 2: (a) Wedgebox plot at SPSla (A bosons only) for 300/fr -1 luminosity, cuts as 
in text. Dotted lines show locations of kinematic endpoints from X234 ~~ > X? decays, (b ) 
Distribution of solutions for SPSla. Here the maximum is at the correct Higgs boson mass 
ofM A ^ 395 GeV. 



Table 1: Relevant masses at the MSSM Test Point I, II (all masses in GeV). 



Particle 


Mass (I) 


Mass(II) 


X? 


119.94 


86.03 


xl 


180.33 


143.09 


xl 


197.98 


166.40 


xl 


317.72 


277.27 


en 


156.17 


127.20 




156.17 


127.20 


A" 


600.0 


700.00 


H° 


609.28 


705.58 



MSSM Test Point I 

fi = 190 GeV M 2 = 280 GeV tan p = 10 M A = 600 GeV 
M^ L R = 150 GeV M^ L R = 250 GeV M u = 1000 GeV. 

which has the mass spectrum shown in Table [U Now the main production modes con- 
tributing to the e + e"/i + /i" signal are the Higgs boson channels H°/A° — » X2X2? X2X3, X2X4; 
'direct' neutralino pair production channels 12 ! X2X3 & X2X4) an d 'mixed' channels 
involving charginos, mainly X2X2 ^ X2 = X2 3 4 production. A random sample of 
e + e~/i + yU~ events^! will therefore not be a clean collection of Higgs boson decays, 

12 As explained in |18j . direct X2X2 channels are suppressed by isospin symmetry, while X3 4X3 4 
are phase-space suppressed. 

13 Note that for MSSM Test Point I the staus have been set more massive than the other 
sleptons to avoid tau-containing decays and generate more decays to the desired leptons. This is 
simply done by hand here for convenience. 



which for the luminosity considered here (300 fbT ) number about 150 against nearly 
800 direct channel and 100 or so mixed channel events. This point then presents a 
more challenging case for applying the mass relation method. 



Test Point I Wedgebox 

200 

180 

160 

140 
S> 120 

(5 

— 100 

'd. 
_k 80 





















%* {•„»*-, ;•" t : 







40 60 80 100 120 140 160 180 200 

M e . e (GeV) 




20 40 60 80 100 120 140 160 180 200 

M« (GeV) 



20 40 60 80 100 120 140 160 180 200 

M« (GeV) 



Figure 3: (a) Wedgebox plot for the MSSM Test Point I defined in the text, for 300 /Z> -1 . 

(b ) origin of each four lepton event labeled as from 'A ', 'H ', direct production 'DP ', or 
otherwise 'Bad '. (c) specific pair of neutralinos decaying to each four lepton event identified 
as xlxl C z % z 4'), X2X3 ('Z2Z3'), x°X° ('Z2Z2'), some other neutralino pair ('Other'), or 
other events 'Bad' 

However, the shape of the wedgebox plot at this point suggests selecting events 
via their decay topology: Fig. [3^ shows a clear 'double wedge' topology implying that 
events due to X2X2 decays are confined to the innermost box bounded by kinematic 
edges at ~ 60 GeV, while those due to X2X3 events are enclosed in the legs of the 
short wedge terminating at ~ 80 GeV. Events from X2X4 decays span both of these 
regions and beyond up to the final kinematic edgd^l at ~ 200 GeV. 

Fig. EIb,c) is color-coded to show the separate distributions of events from differ- 
ent production channels (A, H, or DP 'direct') and by their assorted neutralino-pair 
types; these distributions agree with remarks in the previous paragraph. Note, how- 
ever, the presence of 'bad' events which do not distribute themselves nicely within 
the kinematic bounds and which are typically due to chargino decays. Though nearly 
10 percent of the total number of events are bad events, about half of these are ex- 
cluded by rejecting events outside the overall wedge structure, nicely illustrating the 
strength of the wedgebox technique. 

Without the assistance of the wedgebox plot, one might be led to assume that 
events with M ee MM < 60 GeV correspond mostly to the decay topology of ([T]) with 
i = j = 2. This, however, leads to a Higgs boson mass distribution as shown in 
Fig. UK- There is neither a clear peak nor any kind of distinguishing feature near the 
correct value of Ma = 600 GeV. Evidently, what might seem the natural choice of 
using events from the densest region of the wedgebox plot is not optimal. 

In the MSSM there are other possibilities for these edges aside from those of the form Xj — > 

l + l~Xii including x° — > l + l~X2 an d X2 ~* l + l~xf- Separate studies (including perhaps a Hidden 
Threshold analysis) can probably exclude such possibilities. For parameter set choices examined 
herein, these other potential decay modes have a totally negligible effect. 
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Figure 4: (a) Distribution of solutions for events taken from the inner box region of Fig. 
[3]a; i.e., within the first set of dotted lines at ~ 60 GeV. Here the correct Higgs boson mass 
of Ma ~ 600 GeV does not appear at the peak, (b) Distribution of solutions for events taken 
from the outer legs of Fig. [3]a; i.e., within the regions bounded by dotted lines at ~ 80 GeV 
and ~ 200 Ge V. Here the correct Higgs boson mass of Ma ~ 600 Ge V coincides with the 
peak. 

Instead, events should be selected from the most homogeneous zone of the wedge- 
box plot, which in this case consists of the outermost legs from 80 GeV to 200 GeV, 
corresponding mostly A — > X2X4 (even without looking at Fig. [4b, [18] predicts 
that events in the outer wedge of a double-wedge plot come from Higgs boson de- 
cays). As seen in Fig. Hb, the N-MST method now works splendidly, yielding an 
easily-identified peak at the correct Higgs boson mass, even though now an addi- 
tional sparticle mass plays into the equations. The goodness of fit is more surprising 
considering the number of 'bad' events distributed throughout this zone; inspection 
reveals that these latter, however, often fail to give solutions to equations (jlj)-©, so 
they do not heavily interfere with the shape of the distribution. 

Improvement of the N-MST method via a wedgebox plot, at least for the Higgs 
boson decay topology considered here, is therefore quite straightforward. 

3 The C-MST Method of Cheng et al. 

Next consider the C-MST method. The process treated by Cheng et al. [H] is 
illustrated in Fig. [5] in which the masses m 1; m 2 , m 13 , m 2 4, m 135 , and m 246 must 
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Figure 5: Event topology, taken from [15] 
satisfy a set of equations precisely analogous to ([2])-(j9]): 



m\ 


= pi 


(11) 


m 2 2 


= vl 


(12) 


m n 


= (Pi+Psf 


(13) 




= (P2+P4) 2 


(14) 


™135 


= (P1+P3 + P5) 2 


(15) 


m 246 


= (P2+P4+Pa) 2 


(16) 


P X 

rsum 


= X>f 

i 


(17) 


p y 

r sum 


= 


(18) 



i 



where the transverse momentum sums Psiim are assumed to be calculable from mea- 
surements of associated jet momenta (produced, though not shown, in the gray bub- 
ble in Fig. [5]) and missing momentum (from the LSPs) necessary to balance the 
whole: 

Psum miss Psum miss 

(19) 

1 1 

The specific example considered in [15] was production of two neutralinos via 
squarks in the MSSM, followed by decay via on-shell smuons to muons and two LSPs: 

qq^qq X° X2H #V -> P"VXi) (20) 



giving 



mi — m 2 — m xi 
m 13 = m 24 = m A 

^135 = ™ 24 6 = rn$Q (21) 

In the three-dimensional space of masses (w^o, m^o, m^), each event gives eight equa- 
tions ffTTj) - ffT8]) for the eight unknown LSP momenta p± and p%, assuming the outgoing 
muon momenta P3,Pi,P5,P6 can be measured. The solution to this set of equations 
is again, as for the system of ([2])-®, a quartic equation with 0, 2, or 4 real roots. In 
contrast to the discussion of the last section, however, rather than trying to find the 
point in (m^o, m^o, rap) -space where the density of solutions (= N(m^o, m^o, mpj) 
is maximized, instead the point where the gradient of N is maximized (a heuristic 
argument for this is given in [15]) is sought. 

In [12] this method apparently works quite well at the MSSM parameter points 
studied, giving all relevant sparticle masses to a few percent after collection of data 
samples corresponding to 300 fb _ of integrated LHC luminosity. However, to test the 
robustness of this technique and the extent to which augmentation with the wedgebox 
and/or HT techniques can improve results, consider a new MSSM parameter point 
(see Tab. [T] again for masses): 

MSSM Test Point II 

fi = -150 GeV M 2 = 250 GeV M 1 = 90 GeV 
tan p = 5 Mg^^ = 250 GeV M^ R = 120 GeV 
M A = 700 GeV M~ q = 400 GeV M- g = 500 GeV 



3.1 On- or Off-Shell ? 

First, it should be clear that this method also depends on isolating a fairly pure 
sample of events, in this case with the decay chain fTSOj) . From the wedgebox plot of 
Fig. again generated using ISAJET [23J and the event selection criteria mentioned 
earlier (save no jet cut), one might choose to limit oneself to points M ee wt > 100 GeV 
(away from the "Z-line" at M eew ~ 91 GeV [T7], which is potentially populated 
by Z°Z°-background events as well as undesirable Xj ~~ * Z°Xi events), since the 
resulting subset of the generated events should be mostly of the same XjXj-t°P°l°gy 
(indeed j = 4 hereof]). 

Yet this is not enough; these events must also arise from the decay of on-shell 
sleptons (see Fig. 4 of [T7] for an example of an off-shell wedgebox plot, which is 
visually quite similar), since if the sleptons were off-shell the equalities (TT~3|) and f|T4l) 
would be false. Without the assistance of these equations, the technique of [15] does 
not work at all, since there are not enough constraints to solve for the unknown LSP 
momenta. 



15 Judging from the wedgebox plot alone, this could be interpreted as j — 3 or even a chargino 
decay as previously mentioned; this distinction is irrelevant to the present discussion. 
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Figure 6: Wedgebox at MSSM Test Point II for a luminosity of 90 fb 1 . 

Here the HT method can very easily be applied. As explained in [19], for each 
four-lepton event compute the invariants M 4h M2121, and M« (definitions of 

these reproduced in Appendix A), and plot each of these versus M eejWl (see Fig. [7j). 
Drawing a line at the kinematic edge of M eeiWt ~ 182 GeV, presumed to be equal to 
the mass difference m^o — m^o, the minimum values of each invariant may be read 
off and compared to the expected functions of m^o and m^o (given in Appendix B). 
If this were a genuine off-shell decay, the four sets of solutions would be mutually 
consistent. However, it is easy to check in this case that they are not (the My, 
M31, and Mn minima all imply different negative values of m^o !). Therefore, the 
event sample is consistent with an on-shell decay, and it is safe to apply the C-MST 
methodolog}0. 

3.2 Wedgebox selection 

The wedgebox structural of Fig. [6] is compartmentalized into four substantially 
event-populated regions by the shown red-dashed lines to which will be applied the 

16 A similar analysis can be done if one slepton is on-, one off-shell. 

17 MSSM Test Point II is clearly representative of the general case in the MSSM, where pp 
collisions yield a 'mixed bag' of concomitant neutralino decays. Events on this plot pass the same 
cuts as in Section \TT[ minus the jet cut. 
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Figure 7: Correlations (at 90 fb~ ) among two-lepton and various four-lepton invariant 
masses(see Appendix A for definitions). The off-shell character of decays corresponding to 
the two edges of Fig. at M ee)W ~ 47 GeV and M ee)W ~ 182 GeV may be tested by drawing 
lines at these values and measuring the corresponding minimum values of Mn, M2121, 
and Mu. Since these values do not match what is expected from off-shell formulas the decays 
must instead be on-shell. 



following nomenclature: 

The Zone 1 box, with [M ee and M w ] < 47 GeV, is the most densely-populated region 
of the wedgebox plot and should include all X2X2 events. 

Zone 2 is composed of two rectangles (the legs of wedges) running outwards along 
both axes from the Zone 1 box — satisfying the condition that [M ee and M w ] < 
182 GeV and [M ee orM^but not both] < 47 GeV. Events due to xlxl and xlxl not 
residing in Zone 1 will falfjfl in Zone 2. 

Zone 3, with 47 GeV < [M ee and M MM ] < 182 GeV, lies outside of Zones 1 & 2 and 
should only be populated by X3X3, X3X4 and X4X4 events. 




Figure 8: Mass space scan for MSSM Test Point II, using 1000 random events from the 
wedgebox plot. The correct position of the masses is shown by the dotted lines. 

Consider first using four-lepton events across the whole wedgebox plot — this 
would include several different kinds of neutralino pair events. A scan was performed 
over the (m^o , m^o , m-p) mass space, and the resulting values projected onto the 
(m^o, m^)-plane. This is shown in Fig. [HI from which it is apparent that the gradient 
of N is maximized along the line m^o = m^, which is far from the correct values: the 
C-MST method badly fails in this attempt to find the sparticle masses. Since the 
wedgebox plot clearly indicates that X2X3, X2X4 an d X4X4 production is substantial, 
the failure of the C-MST is not surprising given the inhomogeneity of the data set. 

Now instead we perform an analysis for events in Zone 3 — which should be more 
homogeneous. Yet before proceeding one more feature of the wedgebox plot should 
be taken into account: there is a clear Z-line at [M ee or M w ] ~ 91 GeV. This is due 
to neutralinos decaying to an LSP and a pair of leptons via an on-shell Z° rather than 

18 With the X2X3 (x°Xi) events terminating at M ee or ~ 80 GeV (182 GeV). The endlines 
at ~ 80 GeV are faintly discernible in Fig. however, the forthcoming analysis does not rely upon 
this. 



via a sleptorl^l. The situation is greatly simplified if the events in Zone 3 are further 
curtailed to encompass only those for which 100 GeV < [M ee andM^] < 182 GeV (as 
was done in the previous section). This will exclude the on-shell Z° events along with 
those due to X3X3 or X3Xa- The remaining fairly homogeneous subset of events still 
yields 1000+ signal events (corresponding to 90fb _1 of LHC integrated luminosity), 
80% of which are in fact due to XaXa P arr production (the remainder mostly involve 
colored sparticle decays into the heavier chargino, xf)- With this subset of events 
chosen with the wedgebox plot, a plot analogous to Fig. [8] is generated, as depicted 
in Fig. [9j This time the contours do cluster more closely in the vicinity of the correct 
mass point, though the steepest descent of N is still spread over a line segment of 
the form = m^o + A. 




Figure 9: Same as the previous plot, but now using 1000 events from Zone 3 away from 
the Z-line, i.e., M ee „n > 100 GeV. Contours bunch together closer to the correct point, but 
the point of steepest descent is still ambiguous. 

One reason for persisting inaccuracy is the fact that the momentum conservation 
relation ( fl9l is only approximately true: the MC simulation indicates that each event 
actually has extra px (presumably from undetectable soft parton exchanges) in the 
range —30 GeV > pt ~ 30 GeV (see Fig. [TO] for an example distribution). This is 
formidable in light of the fact that a small inaccuracy in p? (say 2 GeV) can change 
the number of solutions at a given mass point by 10%. Using the correct pt values, 
fits to masses are quite good (easily within 1%). However, it is impossible in practice 
to know the precise pt imbalance, and the best that can be done is to perform a scan 
over this extra pt for each event, taking for instance —10 GeV < px < 10 GeV. This 
is done in making Fig. [TTj, whose maximum does lie somewhat closer to the actual 
LSP mass; however, since the correct value of A = m^ — m^o was assumed in this 1-D 

' ' r Xi 

19 Although the missing energy cut should eliminate most SM Z°Z°, Z°Z°* events, any such 
remnant background surviving would also populate this line. 
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Figure 10: px distribution for one of the transverse directions. Here "real" includes all 
particles (except the LSPs), while "calculated" includes leptons and jets (after the selection 
criteria) only. 

projection, the result is somewhat better than what would be obtained in practice. 
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Figure 11: 1-D projection of Fig. assuming the correct value ofmji — m~o, and scanning 
over —10 GeV < px < 10 GeV. The maximum of the curve, m^o = 98 ± 5 GeV, roughly 
approximates the actual LSP mass (86 GeV), given by the red dashed line. 

The fitted LSP mass in Fig. [11] has a small error but with a slightly up-biased 
central value (primarily due to the fact that the neutralino pairs do not have a fixed 
CM energy). The HT technique can in fact be used to check this result. In [19] it 
was discovered that the four-lepton invariant masses M ih M 2 nu and M 3 ; are 

correlated in a unique way. In particular, a plot of M 3l versus M 4Z x M i2 i or M 3i 



versus M 4; x M 2 m for a suitably large event sample produces a visible curve whose 
shape is strictly controlled by the slepton and LSP masses in the decay chain f[2"Uj) . 
Fig. [12] shows what these plots look like for MSSM Test Point II. Measurement 
of several points along these curves, (i.e., upper envelopes of the shapes formed in 
these scatter plots), yields correct central values for the LSP and slepton masses 
albeit with slightly larger uncertainties (±7 GeV) than those found when just using 
C-MST. 



4 Conclusions 
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Figure 12: HT plots for MSSM Test Point II: points on the top of each curve correspond 
to threshold endpoints which, when matched to formulas in ]1S^ . yield the LSP and slepton 
masses to a few percent. 

Three ways to increase the accuracy and efficacy of the N-MST and C-MST mass 
reconstruction methods have been illustrated: 

1. Use a wedgebox plot to select the most homogeneous sample of events. If the 
events analyzed do not mostly share the same decay topology, both techniques 
fail. If the wedgebox is a simple box, there is nothing to c?o 20 l. However, over 
much of the allowable MSSM parameter space the topology of the wedgebox plot 
is not merely a lone box — if a wedge or composite structure is observed, then 
selecting events from the legs of the wedge or the outer areas proves the most 
effective. Note that this runs counter to the choices made in all the N-MST and 
C-MST publications. In fTSf for instance, some care is taken to describe the 
desirability of so-called 'symmetric events' — where both legs from the original 

20 It is true that at the Snowmass Benchmark point SPSla [20] . a simple box does describe the 
wedgebox topology of a few processes studied at this point. But these SPSla studies [3, [TTl [T3] are 
certainly not representative of many other perfectly allowable parameter set choices and/or signature 
selections. 



parent particle contain the same intermediate particle states. The present work, 
on the other hand, makes the case that the benefits from using un-symmetric 
decay legs, e.g. efficient isolation of events with the same decay chain structure, 
may well trump the convenience of symmetric events in the subsequent MST 
analysis, and therefore unsymmetric events should not be ignored or viewed as 
an unnecessary complication. 

2. Use the Hidden Threshold technique to check decays are indeed on-shell. Since 
the pivotal system of equations or / f77]) -/ f7gj) assume decays through on- 
shell sleptons, it is important to verify this assumption. The HT method can 
be employed to seek a 'null' result to the fit of various J^-lepton invariant mass 
minima. 

3. Scan over the CM in a ±10 GeV-window. Assumptions that the partonic 
CM has no transverse momentum (as implied by equations (TJ|) and / TTgj) ) are 
basically incorrect; while in the N-MST method this does not seem to matter, the 
C-MST method is much more sensitive to this parameter. An 'averaging' over 
Pt improves the result, but perhaps a more detailed analysis should eventually be 
performed as the latest set of structure functions and other knowledge of QCD 
becomes available. 

Note that the time scale required to collect a sufficient number of events to gen- 
erate a wedgebox diagram or apply HT tests is roughly comparable to that needed 
to perform the MST analysis. This is in spite of the fact that the wedgebox and 
HT techniques rely upon populating scatter plots while an MST analysis in prin- 
ciple only requires collection of enough events to simultaneously solve the requisite 
equations. In practice, ambiguity in assigning the leptons and multiple solutions to 
the resulting quartic equation (see bulleted items in Seen. [2] ) as well as experimental 
factors (also enumerated earlier) necessitate a far larger sample of events to perform 
either of the MST analyses discussed. Further, and even more compellingly, without 
augmentation by the wedgebox and HT techniques, applying an MST analysis to 
a quite limited data set is tantamount to wild speculation as to what SUSY chan- 
nels are actually present and the results of such an analysis must be viewed most 
cautiously. 

As noted earlier, the MST analyses presented here assume that the decay chain 
involved is a series of two-body decays via intermediate on-mass-shell sleptons. This 
need not be the case, and this assumption should be tested. The di-lepton distri- 
bution shapes for on- and off-shell decays are not identical [SI E] , and this could 
be used to distinguish between the two possibilities; however, the effects of cuts, 
backgrounds and a finite-sized data set must be considered. Ref. [21] notes that dis- 
tribution shapes for on-shell (sequential two-body decays) and off-shell (three-body 
decays) are effectively indistinguishable for some parameter choices. Also, Ref. [25] 
finds that the shape of the di-lepton distribution may be affected by the nature of 
the neutalinos (the extent to which they are gauginos or higgsinos), illustrating how 
dynamical issues arising from the nature of the coupling involved in a decay may not 
be separable from purely kinematic issues associated with the relevant masses. So 



an alternative to a straight-forward examination of the di-lepton distribution shape 
is desirable [5J- In this study the hidden threshold technique [TS] is used quite 
effectively to fill this need. Herein has been presented (see the C-MST study in Seen. 
[3]) a concrete example of how cross-correlating different invariant mass distributions 
resolves the on-shell vs. off-shell issue (as well as actually improving the accuracy of 
the final masses obtained). Alternative ideas have been put forward in other works. 
In |26p*l a rudimentary sketch of a very Dalitz-esque technique to look for the pres- 
ence of two-body decay chains is presented. This concept is certainly related to the 
HT methodology, and a realistic study applying this idea would be interesting. Refs. 
[9j [27] instead champion a 'Markov chain' approach to analyzing the event sample 
where "no assumption is made about the processes causing the observed endpoints." 
Supposedly then the issue of whether the sleptons are on-mass-shell or off-mass-shell 
is rendered moot to the more modest goal of determining a region of parameter space 
consistent with the data in a non-MST analysis^!. 

The present work may be thought of as an initial installment of a much grander 
programme to fuse all known kinematic mass reconstruction methods together. The 
make-up of this programme consists not only of combining fits from different methods 
for a static event sample set, but also of improving the composition of the event set 
under consideration. In the present work where wedgebox and HT techniques are 
used to select (to the degree possible) events due to a specific x^X°j neutralino pair 
decayingthrough on-shell sleptons^l. Once a fairly homogeneous event sample is 
obtained_j then it becomes straightforward to apply various mass reconstruction 
methods and cross-check them. For example, in the case of Higgs boson decay 
considered here, one could try matching the N-MST results presented here to those 
from a study of 4-lepton invariant mass endpoints [28J - it would be especially 
instructive to compare results at the SPSla parameter point, for example, where both 
techniques, in the total ignorance of sparticle masses, give poor results individually, 
but may give a stronger result in unison [29J. 

An MST analysis is then an attractive option if enough mass-shell conditions 
can be found to match the number of mass components of the invisible final-state 
particles, as is the case for the LSP-generating SUSY decay chains considered herein. 
Further though, the present work shows how endpoint information funneled through 
wedgebox and HT techniques can positively supplement such an MST analysis, as 
in the augmented C-MST study presented in Seen. [31 No mass-reconstruction tech- 
nique is immune from possessing potentially faulty assumptions, and so coupling 
several complementary analysis techniques will in general improve reliability as well 



21 See pages 50-51 therein. 

22 Though suggested, this issue is not explored in any detail in either of these works. 

23 Additional possibilities, e.g. that only one slepton is off-shell, or ^ m^, are similarly testable 
utilizing wedgebox and HT techniques. 

24 An alternative track is attempted in Ref. [27], wherein the idea is to deal with all of the 
complexity of a mixed data set in the mass analysis program, rather than bifurcate the analysis into 
a purifying stage and then an analysis stage. The inherent weakness of this approach is that results 
from studying just the simplest subset of the events are impeded by the need to disentangle more 
confusing events. 



as accurac}@. 

Likewise, consideration of suitable inclusive variables, such as the vari- 
able (30] and its variants [31 J, to augment either an MST study [16J or an end- 
point analysis [TUJ EH [32], [33], [3U [35] has been shown to be beneficial, at least in 
some cases. Then there is the entire array of dynamical (and thus model-dependent) 
information associated with cross-sections and the shapes(spread) of events plotted 
against one(or more) parameters. As noted earlier, kinematics can never really be 
totally divorced from the present dynamics. MSSM/mSUGRA studies combining in- 
formation from cross-sections [27J or distribution shapes [351 [36] with that from an 
endpoint analysis have also been performed and are no doubt the vanguard of many 
more such studies, at least if initial LHC results prove favorable. And, when those 
first major blocks of data from the LHC become available, application of numerous 
analysis techniques — including wedgebox and HT techniques, would be a good idea. 
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Appendix A 

The specific definitions of the 4-lepton kinematic invariants are as follows (lepton 
momenta labeled as in Fig. [T]): 

Ml EE (pi+pi* +P2+P2') 2 (22) 
M 2l2l = {(Pi +PV ~P2 ~P2') A + (Pi +P2> ~P2 -Pl'Y + (Pi +P2 ~ P\' -P2') 4 }/3 
M l2l = {(Pi +VV -V-lf + (Pi +Pl' -P2>) 4 + (Pi +P2 -P2'Y + (P2 +P2> ~ Pl'Y 

+ (Pi - Pi' + P2) 4 + (Pi - Pi' + P2') 4 + (Pi - P2 + P2') 4 + (P2 - P2' + Pi') 4 

+ ipv +P2- Pi) 4 + (Pi' + P2' - Pi) 4 + (pa + P2' - Pi) 4 + (P2' + Pi' - p 2 ) 4 }/12 

M\ x EE {(pi+Pl'+p 2 ) 4 + (pi+Pl'+P2') 4 +(pi+P2+P2') 4 + (p2+P2'+Pl') 4 }/4 

M 4 EE {(p! + p v f + ^ + p 2l f + ^ + p 2 ) 4 + (p 2 + p 2 ,f + (p 2 + pyf + {p V + p 2 ') 4 }/6 
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Appendix B 



The off-shell minima are: 



M 4 f n = ( mj - mi ) 2 + — (23) 



MZ = K - I —P- I (24) 




, 11 + 10^ + 3^ X 

%%n , x | 



M'- = K-mO J. LI (25) 



M,^ = K-mO( ^— ^ | (26) 

For MSSM Test Point II, with rrij - m 1 = 182 GeV, Fig. H yields the following 
values: M™ n = 215 ± 5, GeV, M™ n = 162 ± 5, GeV, M™ in = 170 ± 5, GeV, M™ n = 
115 ± 5, GeV, where generous allowance has been allotted for measurement error. 
Matching each value against the corresponding four formulae above yields, mi < 0, 
m 1 < 102, GeV, mi < 0, and mi < 0, respectively. This is nonsense since the LSP 
cannot have negative mass (moreover, these negative values are all different). 

The decay corresponding to the edge at M eeiW ~ 47 GeV can be similarly ana- 
lyzed, with the same conclusion. 



3 + 2™i + 24 



